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Abstract 

We develop an approach to study the entanglement in two coupled harmonic oscillators. Wc 
start by introducing an unitary transformation to end up with the solutions of the energy spec- 
trum. These are used to construct the corresponding coherent states through the standard way. To 
evaluate the degree of the entanglement between the obtained states, we calculate the purity func- 
tion in terms of the coherent and number states, separately. The result is yielded two parameters 
dependance of the purity, which can be controlled easily. Interesting results are derived by fixing 
the mixing angle of such transformation as £ . We compare our results with already published work 
and point out the relevance of these findings to a systematic formulation of the entanglement effect 
in two coupled harmonic oscillators. 
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1 Introduction 



Entanglement is one of the most remarkable features of quantum mechanics that does not have any clas- 
sical counterpart. It is a notion which has been initially introduced and coined by Schrodinger [1] when 
quantum mechanics was still in its early stage of development. Its status has evolved throughout the 
decades and has been subjected to significant changes. Traditionally, entanglement has been related 
to the most quantum mechanical exotic concepts such as Schrodinger cat [I], Einstein- Podolsky- Rosen 
paradox [2] and violation of Bell's inequalities [3]. Despite its conventional significance, entanglement 
has gained, in the last decades, a renewed interest mainly because of the development of the quantum 
information science [3]. It has been revealed that it lies at the heart of various communication and 
computational tasks that cannot be implemented classically. It is believed that entanglement is the 
main ingredient of the quantum speed-up in quantum computation [3]. Moreover, several quantum 
protocols such as teleportation, quantum dense coding, and so on [5j [6j [7J O [TOj, [11] are exclusively 
realized with the help of entangled states. With this respect, many interesting works appeared dealing 
with the development of a quantitative theory of entanglement and the definition of its basic measure. 
These concern the concurrence, entanglement of formation and linear entropy |12t 1131 HH1 115j. 

Entangled quantum systems can exhibit correlations that cannot be explained on the basis of 
classical laws and the entanglement in a collection of states is clearly a signature of non-classicality 
|16| . Furthermore, in the last few years it has become evident that quantum information may lead to 
further insights into other areas of physics [T7J. This has led to a cross-fertilizing between different 
areas of physics. It is worthy of note that the nonlinear Kerr effect [18] has been considered as the 
most famous source of physical realization of photon pairs of entangled polarization states. However, 
it raises a number of difficulties to the control of photons that are traveling at the speed of light. This 
is why so much attention has been paid recently to the entangled states of massive particles as they 
are viewed to be much more easy to control [I7J Q2] . 

On the other hand, the harmonic oscillator machinery plays a crucial role in many areas of physics. 
These are the Lee model in quantum field theory [20] , the Bogoliubov transformation in superconduc- 
tivity [21] . two- mode squeezed states of light [22 ^ I23 t 12^]. the covariant harmonic oscillator model for 
the parton picture [25], and models in molecular physics |26j. There arc also models in which one of 
the variables is not observed, including thermo-field dynamics [27], two-mode squeezed states [2"8"l [2T?] . 
the hadronic temperature [30], and the Barnet-Phoenix version of information theory [31] . These 
physical models are the examples of Feynman's rest of the universe. In the case of two coupled har- 
monic oscillator, the first one is the universe and the second one is the rest of universe. For sake 
of the mathematical simplicity, the mixing angle (rotation of the coordinate system), in the above 
mentioned references, is taken to be equal ^. This means that the system consists of two identical 
oscillator coupled together by a potential term. 

In the context of the entangled massive particles, we cite the recently achieved investigation of a 
specific realization of two coupled harmonic oscillator model by the authors of reference [T^]. In fact, 
they calculated the interatomic entanglement for Gaussian and non-Gaussian pure states by using 
the purity function of the reduced density matrix. This allowed them to treat the cases of free and 
trapped molecules and hetero- and homonuclear molecules. Finally, they concluded that when the 
trap frequency and the molecular frequency are very different, and when the atomic masses are equal, 
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the atoms are highly-entangled for molecular coherent states and number states. Surprisingly, while 
the interatomic entanglement can be quite large even for molecular coherent states, the covariance 
of atomic position and momentum observables can be entirely explained by a classical model with 
appropriately chosen statistical uncertainty. 

Motivated by the mentioned references above and in particular |19| . we undertake to develop a 
new approach to study the entanglement in two coupled harmonic oscillators. It is based on a suitable 
transformation having the merit of reducing the relevant physical parameters into two: the coupling 
parameter 77 and mixing angle 9. It turns out that we can easily derive the solutions corresponding to 
the energy spectrum. Then, the obtained solutions are used to construct the coherent states through 
the standard method. In order to characterize the degree of entanglement, we calculate, within the 
framework of the coherent states, the purity function. Then the final form of the purity is cast in 
terms of 77 and 9. Our finding shows two interesting results: the first one tells us that the present 
system is not entangled at 77 = 0, as expected, and highly entangled at large 77 (Figure 1). The second 
one is when we fix 9 = ^, the purity behaves like the inverse of cosh 77 and the corresponding plot 
(Figure 2) shows that the purity is ranging between and 1. It is worthy of notice that, in this case, 
the purity becomes one parameter dependent, which means that it is easy to control. 

Subsequently, we evaluate the purity in terms of the number states. In doing so, we use the well- 
known relation to express the number states \ni,ri2) as function of the corresponding coherent states 

1 a, /?). Then after a lengthy but a straightforward algebra, we end up with the final form of purity. To 
be much more concrete, we restrict ourselves to some interesting cases that are (n\ = l,n2 = 0) and 
(ni = l,r&2 = 1). For the first configuration, the obtained purity is simply a ratio of hyperbolic and 
sinusoidal functions, which tells us that the entanglement is maximal at large 77 for all 9 (Figure 3). In 
the particular case when 9 = ^, the purity is typically a ratio of a hyperbolic cosine function, which 
shows clearly that the purity is positive as it should be (Figure 4). The second configuration gives also 
a mixing dependence between the hyperbolic and sinusoidal functions where the corresponding plots 
(Figures 5 and 6) show some difference in the form with respect to the first one. In both cases, we 
notice that the numerators are always hyperbolic cosine of even 77 and denominators are also power of 
the function cosh 77. 

The present paper is organized as follows. In section 2, we review the derivation of the solutions 
of the energy spectrum of two coupled harmonic oscillators |32] , These will be used to build the 
corresponding coherent states and therefore evaluate the purity function of the reduced matrix elements 
in section 3. The final form of purity function is subjected to different investigations where we underline 
its dependence to two physical parameters 77 and 9. In section 4, we evaluate the purity in terms of 
the number states after a series of transformation. Two interesting case of the purity will be discussed 
in section 5. Finally, we give conclusion and perspective of our work. 

2 Energy spectrum solutions 

In doing our task, we consider a system of two coupled harmonic oscillators parameterized by the 
planar coordinates (Xi,X2) and masses (mi, 777,2). Accordingly, the corresponding Hamiltonian is 
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written as the sum of free and interacting parts 

Hi = ir— p i + 7T— p 2 + \ {C x Xl + C 2 X\ + C 3 X 1 X 2 ) (1) 
lm\ 2m 2 2 ' 

where Ci,C 2 and C3 are constant parameters. After rescaling the position variables 

x\ = fj,X\, x 2 = H~ 1 X 2 (2) 

as well as the momenta 

Pl=fJ,~ 1 Pl, P2 = ^P2 (3) 

Hi can be written as 

H 2 = — {pi + pi) + - (cix\ + c 2 x\ + C3X1X2) (4) 
where the parameters are given by 

ji = (mi/m 2 ) 1/4 , m = (mim 2 ) 1/2 , ci = C\J — , c 2 = C 2 .—, c 3 = C 3 . (5) 

V Tni y m 2 

As the Hamiltonian (j3D involves an interacting term, a straightforward investigation of the basic 
features of the system is not easy. Nevertheless, we can simplify this situation by a transformation to 
new phase space variables 

Va = M ab x b , p a = M ab p b (6) 

where the matrix 

r** \ I cos l -sin I \ 

V Sm 2 COS 2 / 

is a unitary rotation with the mixing angle 9. Inserting the mapping © into one realizes that 6 
should satisfy the condition 

tanfl= ° 3 (8) 
C2 - ci 



to get a factorizing Hamiltonian 



H 3 = ^(pi+P 2 2 ) + t(e%l + e-^ 2 ) (9) 



where we have introduced two parameters 



k = Jcic 2 -4/4, & = £l±WEE g±g (1Q) 



2k 

under the reserve that the condition 4cic 2 > c| must be fulfilled. The parameter 77 is actually measuring 
the strength of the coupling. 

For later use, it is convenient to separate the Hamiltonian ([9]) into two commuting parts and then 
write H3 as 

H 3 = e^Ui + e^U 2 (11) 

where 7~Li and H 2 are given by 

Hi = ±e-"ft + feV, U 2 = ^p\ + l -e^yl (12) 
2m 2 2m 2 
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One can see that the decoupled Hamiltonian 

Ho = — Pi + -y? + — pi + -yl (13) 
u 2m yi 2 yi 2m y2 2 y2 v ; 

is obtained for 77 = 0, which is equivalent to set C3 = 0. 

The Hamiltonian H3 can simply be diagonalized by defining a set of annihilation and creation 

operators. These are 



2hw 1 V2fc 



-e 2 j/i + — ; e zpj, = W— j/j e 2 (14) 



with the frequency 



V m 

and e = ±1 for i = 1,2, respectively. They satisfy the commutation relations 

[aj,at] = <% (16) 
whereas other commutators vanish. Now we can map H3 in terms of a, and aj as 

i?3 = faj ^e^ajai + e~ T, a» l a2 + cosh 77^ . (17) 
To obtain the eigenstates and the eigenvalues, one solves the eigenequation 

H 3 \n 1 ,n 2 ) = £ ni ,n 2 \ni,n 2 ) (18) 

getting the states 

(a\) ni (ai) n2 , . , x 

ni,n 2 ) = ^4=^10,0) (19) 

as well as the energy spectrum 

-E , 3,ni,n 2 = ^ {e r 'ni + e"' ? n2 + cosh 77) . (20) 

It is clear that these eigenvalues reduce to those of the decoupled harmonic oscillators, namely 
Huj (ni + n 2 + 1). This shows clearly that the presence of the coupling parameter 77 will make dif- 
ference and allow us to derive interesting results in the forthcoming analysis. 

To show the correlation between variables, let us just focus on the ground state and write the 
corresponding wavefunction in y-representation. This is 



My) - <yi,y 2 |o,o) = y^exp{-^ (e^ + e^y 2 )} (21) 

which can easily be used to deduce the ground state wavefunction in terms of the variables (xi,x 2 )- 
Therefore, from the unitary representation we find 

tpo(x) = (xi,x 2 \0,0) 

rru.> r / e 9\ 2 _ / e #\ 2 1 1 . . 

xi cos - - x 2 sin - I + e 'I x\ sin - + x 2 cos - J > . (22) 

We notice that f)21|) is separable in terms of the variables y\ and y 2 , which is not the case for (|22p in 
terms of x\ and x 2 . We close this part by claiming that the obtained results so far will be used to 
study the entanglement in the present system. 







V^T exp< 


~~2h 
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3 Entanglement in coherent states 



As we claimed above, we implement our approach to study the entanglement of two coupled harmonic 
oscillators. Actually, it can be seen as another alternative method to recover the results obtained 
in [19] not only in a simpler way but also with less physical parameters of control. To start let us first 
introduce the coherent states corresponding to the eigenstates |rai,n 2 ) given in (|19|) . As usual, we can 
use the displacement operator to define the coherent states in terms of two complex numbers a and 
p. These are 

\a, p) = D(a 1 ,a)D(a 2 ,P)\0, 0) (23) 



which gives the wavefunction 



^afi (2/1,2/2) 



A1A2 



7T 



1/2 



exp 



A 2 



a 



a 



+ y/2a\i 



2/1 



where we have set the quantities 

Ai = 



2 



rrik^ 



A 2 



_2 

e 2 



a 2 

o 2/2 



+ V2pA 



22/2 



(24) 



(25) 



In terms of the original variables (Xi,X 2 ), (|24p reads as 

1/2 



$ Q/3 (X 1 ,X 2 ) 



7T 

x exp 
x exp 



exp 



4 LucosfXi - isin|X 2 (^sinfXi + icoslXs 



V2a\i (fi cos |Xi - 1 sin f X 2 ) + \/2/3A 2 ( // sin fXi + ± cos §X 2 
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\§l 
2 



if 
2 



(26) 



As it is clearly shown in the wavefunction (|26p , the non-separability of the variables will play in crucial 
role in discussing the entanglement in the present system. This statement will be clarified later on 
when we will come to the analysis of the role of the involved parameters. 

At this level we have set all ingredients to study the entanglement in the present system. All we 
need is to determine explicitly the purity function that is a trace of the density square corresponding 
to the obtained eigenstates. More precisely, we have 



P = TV 

which in terms of the above coherent states reads as 



a/3 



dX x dX' x dX 2 dX' 2 $> afi (X 1 ,X 2 ) §* afi {x[,x 2 ) <D Q/3 {X[,X 2 ) $* a p (Xt, x' 2 ) 



(27) 



(28) 



Upon substitution, we obtain the form 



A1A5 



7T 



dX 1 dX[dX 2 dX' 2 e 



V(A? cos 2 §+A 2 sin 2 |)(x 2 +X; 2 )-4(A 2 sin 2 f+A 2 cos 2 § )(Xf+X^ 2 ) 



1 

xe 2 



xe 



(A 2 -A 2 ) S ine(X[X 2 +X 1 X 2 +X[X! 2 +X 1 X! 2 ) ^(^Ai cos § + ^A 2 sin §) (X 1+ X[ ) 



(29) 
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This integral can easily be evaluated by introducing an appropriate transformation. This can be done 
by making use of the following change of variables 





Xy 








X'l 




1 




x 2 




~ 2 


V 


x' 2 


) 





V 



0)1 


A 4 1 
-A* 


0)1 


fj,\/l— 2a 

Wl 


/Vl+2a 

0)1 


fJ-\/l— 2a 




flUJ2 





^0)2 


y/l-2a 


Vl+2a 


flUl-2 





£taj 2 


y/l-2a 


Vl+2a 






-y/2fj,uj 2 



\ 





( 


»1 


\ 






u 2 














\ 


U4 


J 



(30) 



where u)\, d 2 and a are given by 
1 



Co'! 



A 2 cos 2 § + A 2 , sin 2 f 



W2 



A 2 sin 2 | + A 2 cos 2 | 



a = — - (A 2 — A|) sin^wia;2. (31) 



By showing that the determinant of such transformation is , it is easy to map P a p in terms of the 
new variables as 



a/3 



7T- 



;XlX 2 UlU 2 e 



-2|a| 2 -2|/3| 2 -a 2 -a* 2 -/3 2 -/3* 2 



+oo 



du\ e 



duz e 



u i+ 77\ 2a [ Al ( Wl cos l +W2 sin |)(a+a*)+A 2 (o)i sin |-o> 2 cos |)(/3+/3*)]«i 



m§+ /jv 2a [Ai(o>i cos § -o) 2 sin |)(a+a* )+A 2 (wi sin |+w 2 cos |)( / 3+^*)]u 3 



(32) 



Performing the integration to end up with the result 

Pa/3 M) = 



2 cosh 277 sin 2 | cos 2 | + cos 4 | + sin 4 



(33) 



This is among the interesting results derived so far in the present work. Indeed, it shows clearly that 
the purity depends on the physical parameters (77, 6) rather than the complex displacements (a, f3) 
and hereafter it will be denoted by P(r/,9). Furthermore, the obtained purity is two parameters 
dependent, which means that it can be controlled easily. If one requires the decoupling case (77 = 0), 
P(rj, 9) reduces to one as expected and therefore there is no entanglement. 

To understand better the above results, we recall that the purity is related to linear entropy by 
the simple form 

L = 1 — P (34) 

where P lies in the interval [0,1]. Now let us proceed to plot the purity for a range of 77 and by 
considering 9 € [0, ir). From Figure 1, it is clear that the purity, as function of 77, is symmetric with 
respect to the decoupling case rj = 0. It is maximal for 77 = 0, which really shows that the system 
is disentangled. After that it decreases rapidly to reach zero and indicates that the entanglement is 
maximal. More importantly, the purity becomes constant whenever 9 takes the value zero or ir. This 
behavior of the purity traced in below tell us that one can easily play with two parameters to control 
the degree of entanglement in the present system. 
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Figure 1: Purity in terms of the coupling parameter rj and the mixing angle 9. 

Specifically at 9 = ^, we obtain a simple form 

/ 7T 



1 



(35) 

cosnr/ 

which is one parameter dependent and can be adjusted only by varying the coupling ij to control the 
degree of the entanglement. To be much more accurate, we underline such behavior by plotting ([35]) 
in Figure 2: 

P(77) 





-10 -5 5 10 

Figure 2: Purity in terms of the coupling parameter r] for the mixing angle 9 = | . 

From the above figure, one can deduce two interesting conclusions. The first one tells as that P(r], 9) 
is bounded, i.e. < P < 1, as expected. The second one shows clearly that the purity goes to zero 
for a strong coupling, which indicates the entanglement is maximal. 



4 Entanglement in number states 



To gain more information about the behavior of the present system, we evaluate the degree of the 
entanglement between inter states. For this, we consider the relation inverse to express the number 
states in terms of the coherent states. This is 



ni,n 2 ) 



e H 2 e l< 



a, P) 



(36) 



a=0,/3=0 
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In the y-representation, (|36p leads to the wavefunction 

(919 9 19 

/xcos-Xi - -sin-X 2 ,/isin-Xi + - cos -X 2 ) = $ ni n 2 (Xi,X 2 ) 



1 d ni d n2 J£p , v v , 

-e 2 e 2 ^ a p{X u X 2t 



0nln 2 ! <9a n i dp 



11-2 



(37) 



a=0,/3=0 



where Q a p (Xi,X 2 ) is given in (|26p . This will be implemented to study the purity in terms of the 
number states and discuss different issues. 

Returning back to the purity definition, we have 

P nin2 = I dX x dX' x dX 2 dX' 2 $> nxn2 (X U X 2 ) $; in2 {X[,X 2 ) $ nm2 (X[,X' 2 ) K in2 {X U X' 2 ) . (38) 



Using (|37p to obtain the form 



dXidX[dX 2 dX 2 



ni\n 2 \ 



g n i g n 2 ^il 



s™i a"2 Jj^Jlst 



1 



I 



e 2 e 2 $ Ql(3l Uicos-Xi sin -X 2 , fj, sin -X x + - cos -X 2 



/i 2 



8«i a"2 ^l^aJi^iMf 



1 



e 2 e 2 $* „ lyucos-X^ sin -X 2 ,/j sin -Xj H — cos -X 



/' 



Q"2 "4l 2 l^4| 2 



1 



e 2 e 2 $ Q3/ 3 3 //cos-X| sin -x' 2) /xsin -X{ + - cos -x 2 



X 8a? 



e 2 e 2 $ 



ju 2 
1 



oti,/3i=0 

(39) 

«2^2*=0 



"3,^3=0 



\i 2 

After some algebra, we show that the purity takes the form 



cwjfc I ^cos-j;i - -sm-i 2 ^sm-Ii + - cos -X 2 j 



Pn 



( AlA2 ) 2 I dXldX[dX2dx ' 2 e-^^^^m 



\2 

21 

2 



\TTmWJ tJida^dff 2 

(/j, cos §Xi - 1 sin | X 2 ) 2 + (fJ, cos §X( - A sin §X£ V + (/i cos §X( - A sin § X 2 ) + (u cos f Xi - A sin §X£ Y 
(,usin§ Xi+i cos f Jf 2 ) 2 + (/i sin§X(+A cos f X^) V (/i sin §Xi+A cos f X/,) 2 + (/j sin §X(+A cos §X 2 Y 



\2 

23. 

2 



? V2/ i [(A 1 (a 1 +a4)cos|+A2(/3i+/34)sm|)Xi + (A 1 (a2+a3)cos|+A 2 (/3 2 +/3 3 )sinf)X(] 
[(A 1 (a 1 +a 2 )sin|-A 2 (/3 2 +/3 1 )cos|)X2 + (A 1 (a3+Q4)sinf-A 2 (/33+/3 4 )cosf)X^ 



e f 



(40) 



This can be written, in a compact form, as 

AiA 2 



nin 2 



2 4 <9 ni 9™ 2 



n 



where z T = I X\ X[ X 2 X' 2 ) , the matrix ^4 is given by 



<FZe 



-Z t -A-Z+B t -Z+C 



A 



( 


A u 





A 13 


A13 









An 


A 13 


A13 








A13 


A 33 







{ 


A 13 


A13 





^33 





(41) 



(42) 
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such that their components read as 



An = H ( K cos z - + \\ sin z - 



A 



33 



1 

/' 



\2 ■ 2 ° i \2 2 17 

A x sin - + A 2 cos - 



.4 



13 



^-^sintf (43) 



and the matrix B takes the form 

^ /xAi (ai + 04) cos I + /^A 2 + Pa) sin | ^ 
AiAi (a 2 + a 3 ) cos | + fi\ 2 (/3 2 + /3s) sin § 
^A 2 + /3 2 ) cos § - iAi (ai + a 2 ) sin | 
V ^A 2 (/3 3 + /3 4 ) cos § - ^Ai (a 3 + a 4 ) sin f ^ 



B = V2 



(44) 



To go further in evaluating the purity, we perform a method to simplify our calculation. This can 
be done by introducing the change of variables 



/ 


Xi 








X'i 








x 2 






V 


x' 2 


) 





I 









2^\~2a 










2^\-2a 


2^t^l+2a 


tiu> 2 







2Vl-2a 


2y/l+2a 


/J,u>2 







2Vl-2a 


2^1+2a 






^fJ,UJ 2 
V2„, _ 





( 


XI 








X2 








X3 






\ 


X4 


) 



(45) 



where the corresponding measure is dXidX^dX^dX^ = Jdxidxidx 2 dx 2 and the Jacobian reads as 



J= / 

AiA 2 ^(A? cos 2 § + A 2 , sin 2 f ) (A? sin 2 § + A 2 cos 2 §] 
This performance allows us to map (j41[) as 



P, 



niri2 



AiA 2 
irni\n 2 \ 



4 ^ni gn2 



J»2 



where Q* 



I? = a/2 



xi x 2 x 3 x 4 I and D is the transpose of D, such as 



/ ui cos |+w 2 sin § . , , \ , wisin|-aj 2 cos| . ,„ Q . a , o\ 

1 - -Ai (ai + a 2 + a 3 + a 4 ) H \_ 0n — -A2 (p\ + /3 2 + /3 3 + P4J 



^wiAi (ai + a 4 - a 2 - a 3 ) cos § + ^-ui\ 2 (p\ + /3 4 — /3 2 — /3 3 ) sin ■ 



2 y-^l I »-t4 *-*o y 2 I 2 w^/^ 1 ^4 pz y am 2 

wi cos |~wa sin (|) \ / , , 1 \ 1 sin |+W2 cos I . ,„ „ „ , o \ 

^Ai (ai + a 4 + a 2 + a 3 ) H 2 , — ^-A 2 (p\ + /3 4 + /3 2 + /3 3 ) 



V 



2^1+2^ 1-41^1 -dy 1 2^1+2^ 

-^w 2 Ai (ai + a 2 - a 3 - a 4 ) sin | + ^u 2 X 2 (/3i + /3 2 - /3 3 - /3 4 ) cos | 



2 "/z^vi ^uii 1 pm 2 I 2 \f^l I A^Z A^d ^4; ^^'J 2 

Since the above integral is Gaussian, then after some algebra we end with the form 



P, 



nin 2 



AjA 2 
irni\n 2 \ 
2t 



4 Qni Qri2 

J BiWW exp 



(46) 



(47) 



(48) 



2s 



u 2 2w 2u 2w 2s 

-ai H aiat 2 aia 3 H aia 4 H aipi 

P P P P P 

2w 2u 2t „ 2s 



Lb U 

«i/3 2 -ai/33 H ai/3 4 H — a 2 H — ~ a 2«3 ^«2«4 «2^i H — -ot 2 (3 2 

P P P P P P P P 

2t 2s „ u 2 2v 2s „ 2t . 2s . 2t 

H — a 2 pz a 2 p4 + —a 3 H a 3 a 4 a 3 pi H a 3 p 2 H a 3 p 3 a 3 p 4 

P P P P P P P P 

u 2 2i 2s 2t 2s ti o 2tt) 
H — a 4 H a 4 pi a 4 /5 2 a 4 /5 3 H a 4( o 4 & H Pip 2 

P P P P P P P 

2u n _ 2v u o 2v 2w u ^2 2w u . 

P P P P P P P P _ 



(49) 
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where we have set the involved parameters as 

, mk 

p = V 

2mk 



2 cosh(2r/) + cot 2 - + tan 2 - 



2mk 



sinh 2r] 



h 2 

mk cosh rj 



cosh(2r?) + tan 



ir 



mk 
-4—5" sinh 77 
nr sm 



2mk 
= ~W 
cos 6 



n 2 

cosh(2r/) + cot 2 



(50) 



h 2 s'm8 

We are still looking for the final form of the purity, which can be obtained by calculating the partial 
derivatives. These can be performed in different ways and may be it is easier to proceed step by step. 
Indeed, we factorize the exponential function and then map each factor into a series expansion. This 
operation has been postponed to Appendix A and the yielded result is 



P, 



nin 2 



M) 



I) 2 Vi!«2<) 2 



E 



sin(0) y2 cosh 2j?+tan 2 ( § ) +cot 2 ( § ) i+j+k+l+r=2 ( ni +„ 2 ) 

where the coefficients C nin2 are given by 



a 



nin 2 



(i,j,k,l,r)u i v j w k t l s r (51) 




1 



2~«4 I3+U— is+'e— J7-H - — »"l-|-r3— »'5+* , 6— »"7+i2— Cl— C2 



7 r e -i 

n e 

e=lr ( 



(r e _i - r e )! 



(52) 



(i - (ii - i 2 )! («2 - ^3)! (*3 - i4)^7^7!c3!c4!c 5 !c6!c7!c 8 !c 9 !cio! ' 
It is clear that the final form of the purity is actually only depending on two parameters, i.e. 77 and 
6. On the other hand, it is easy to check that P ni „ 2 is symmetric under the change of the quantum 
numbers n\ and 



5 Two special cases 

To be much more accurate let is illustrate some particular cases. With these we will be able to get 
more information from the above purity about the degree of entanglement. In the beginning, let us 
choose the configuration (m = 0, n 2 = 1), which means that we are considering now the entanglement 
between the ground state of the first oscillator and the first excited state of the second one. In this 
case, (I5T1) reduces to the form 



PoiM) = ^ . — E C i(i,j,k,l,r)uVw k t l s r (53) 

sin 9^2 cosh(2r;)+tan 2 j+cot 2 5 i +r+ j + f.+i=2 

which can be evaluated to obtain 

P01 (V, 0) = , KpJ (u 2 + v 2 + w 2 ) (54) 

sin6»W2cosh(2r))+tan 2 |+cot 2 § 

and after replacing different parameters, one gets the final result 

, 3cosh(4r7)+4(tan 2 | + cot 2 |)cosh(2r/) + 2tan 4 | + 2cot 4 | + l , 

P01 (V, 0) = — ^ 2 2j — ^ 2 ■ (55) 

sin 9 (2 cosh(2r/) + tan 2 | + cot 2 f ) 2 

This is a nice form, which can be worked more since it is only function of two physical parameters r\ 
and 0. Indeed, we plot it in Figure 3: 
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Figure 3: Purity Poi as function of the coupling parameter rj and mixing angle 9 for the quantum 

numbers (n\ = 0, n,2 = 1). 

Here we have the same conclusion as in Figure 1 except that the present plot is showing some de- 
formation at the point r\ = 0. Otherwise, for certain values of 9 the purity is not always holding a 
maximum value at 7] = 0. More precisely, at this point it decreases to reach 1/2 at 9 = \ and then 
increases to attends 1 at 9 = ir. This is because in the present case the masses are equal and the same 
conclusion is obtained in [19]. 

Now let us look at some interesting situations by fixing the mixing angle 9 and varying the coupling 
parameter rj. In particular when 9 = ^, Poi reduces to the form 

7T \ _ 3 cosh(47?) + 8 cosh(27?) + 5 
2 ) ~ ~ 32 cosh 5 7] ~' 

This can be plotted to obtain Figure 4: 

P(77) 



(56) 




77 

-10 -5 5 10 

Figure 4: Purity P01 as function of 77 measuring the entanglement between the ground state n\ = and 

the first excited state 712 = 1 for 9 = |. 

Compared to Figure 2, we notice that the behavior of the purity in terms of the coupling parameter r\ 
is large. As long as r/ is large the entanglement is going to hold the maximum value. It shows clearly 
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the role playing by 77 and thus allows an easy control of the degree of the entanglement. This may 
give some hint about an experiment realization of the present case. 

Now let us look at the case of the entanglement between the two first excited states of the two 
oscillators, i.e. n\ = ni = 1. This result gives 



Pn (v,o) 



2 4 



sin Q\Jl cosh(2»;)+tan 2 |+cot 2 | j_|_j_|_/ c _|_;_|_ r= 4 

after lengthy but simple calculations, we find 



£ C u (i,j,k,l,r)uVw k t l s r 



(57) 



Pi 



ii 



sin 0^/2 cosh(2r;)+tan 2 |+cot 2 § 



(58) 



x (u 4 + v 4 + w 4 + 2 2 s 4 + 2 2 t 4 + 2 4 sH 2 + 2u 2 v 2 + 2u 2 w 2 + 2v 2 w 2 + 2 4 ustv - 2 4 ustw 

n4 2 2 n4 2j.2 o4j.2 2 n4 2 2 n4j.2 2 o4 2 2 , o3 2 , n3 j.2\ 

—2 u s — 2 u t — 2 t w — 2 s w — 2 t v — 2 s v + 2 vws + 2 vwt j . 
Finally, we obtain 



4 sin 6 [2cosh(2r/) + tan 2 f + cot 2 f] 2 



9 cosh (877) + 16 ( tan" — + cot" — ) cosh 677 



( 96 tan 4 - + 96 cot 4 - - 36 j cosh (4r?) + 240 f tan 2 - + cot 2 ) cosh (2??) 



+8 tan 8 - + 8 cot 8 - - 64 tan 4 - - 64 cot 4 - + 459 
2 2 2 2 



(59) 



Comparing this with (J55J), we notice that the numerator of both of them is containing a hyperbolic 
cosine function of a even number of coupling parameter r\ and the denominators are power of cosh 77. 
To go further, we plot (I58|) in Figure 5: 



P (rj, 9) 




Figure 5: Purity Pn as function of the coupling parameter 77 and mixing angle for the quantum 

numbers (ri\ = 1,77,2 = 1)- 

Clearly, we see that for certain values of 9 the purity is not always holding a maximum value at 
decoupling case, i.e 77 = 0. At this point, the purity decreases to reach 1/2 at 6 = ^ and then increases 
to attends 1 at 9 = n. 
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Furthermore, (|59p can be worked much more to underline its behavior. The simplest way to do so 
is to fix the mixing angle 9 and play with the coupling parameter rj. For instance, by requiring 9 = ? 
we end up with the form 



Pn(ri,0) 



9cosh(8r/) + 32cosh(6r/) + 156cosh(4r/) + 480cosh(2r?) + 347 



(60) 



2048 cosh 9 r] 

This shows clearly that P\\ {rj, 9) is one parameter dependent and therefore it can be manipulated 
easily. For more precision, we plot (|60p in Figure 6: 

P(77) 




-10 -5 5 10 

Figure 6: Purity P\\ as function of rj measuring the entanglement between the first exited states 

(m = l,n 2 = 1) for 9 = f . 

This is showing a difference with respect to Figure 4. It is clear that as long as rj is small the purity 
increases rapidly to reach its maximal value. Also it decreases rapidly to attend zero for large 77, which 
means that the system is strongly entangled. 



6 Conclusion 

The present work is devoted to study the entanglement of two coupled harmonic oscillators by adopt- 
ing a new approach. For this, a Hamiltonian describing the system is considered and an unitary 
transformation is introduced. With this latter, the corresponding solutions of the energy spectrum 
are obtained in terms of the coupling parameter r\ and the mixing angle 9. It is clearly seen that when 
r\ = 0, the system becomes decoupled and therefore nothing new except harmonic oscillator in two 
dimensions. 

To study the entanglement of the present system, we have introduced the purity function to 
evaluate its degree. In the beginning, we have realized the corresponding coherent states by using the 
standard method based on the displacement operator. These are used to determine explicitly the form 
of the purity in terms of the physical parameters rj and 9. Also, the obtained result confirmed the 
range of the purity that is < P < 1. Moreover, we have clearly shown that purity is easy to control 
and can also be cast in a simple form when we fix 9 = ^. In such case the purity is obtained as the 
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inverse of the hyperbolic function cosh 77 and the disentanglement simply corresponds to switching off 
V- 

Subsequently, we have used the relation inverse between the number of states and the coherent 
states to determine the purity. After making different changes of variable, we have got a tractable 
Gaussian form, which was integrated easily. The final result showed that the purity is two parameters 
dependent. This allowed us to illustrate our finding by restricting ourselves to two particular cases. 
In the first configuration, we have considered the entanglement between the ground state and excited 
state, i.e. (m = 0, 77.2 = 1) where the purity is exactly obtained. In the second configuration we 
studied the entanglement between the states (n\ = 1,77-2 = 1). In both cases, we have analyzed the 
case where = 5, which showed that a strong dependence of the purity to the hyperbolic cosine 
function of even coupling parameter. 

On the other hand, the system of two coupled oscillators can serve as an analog computer for many 
of the physical theories and models. Therefore, one can extend the method developed here to study 
the entanglement in other interesting systems those illustrating the Feynman's rest. Furthermore, one 
immediate extension is to consider the case of a coupled systems submitted to an external magnetic 
field. This work and related matter are actually under consideration. 
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Appendix A: Final form of purity 



In this appendix, we show how to derive the final form of the purity given in ()5ip . Indeed from ([49 P . 
we obtain the result 



AiA 



1^2 



7T77i!772! 




J! 



J E 

i,j,k,l,r=0 \e=li e =(K''e-l 
1 r e -i 

E 



E 



1 



E 



1 



j e =0Ue-l _ Je)! fc e =o(^e-l 



V 

l e =G\le-l ~ l e )!r e =o(>"e-l ~ r e )\ 



d n i ai 
-a 



(9"2 a a 7 \ ( d n 2 a a 8 



(a i ,/3 i ) = (0,0) 




fee)! 
2 i-i 4 



2«y (2w_y (2ty (2£ 

PJ \ PJ \PJ \pJ hiW-fo^kW- 



da^ 1 4 



9"2 qa 5 \ ( d n 2 ga 6 \ 

e^ 2 Pl J Va/32 2P2 J 



(Al) 



14 



where different parameters are given by 



ai = 2in + h — ii + h — l 2 + l — h + r\—r 2 + r — r\+j 3 + k — fa 

a 2 = 2ii - 2iii + 12 ~ h + h + k ~ h + r 6 - r 7 + r 5 - r 6 + j 3 + k 3 

a 3 = 2i 9 - 2ii + «3 - «4 + ^5 - + h - h + r 7 + r 4 - r 5 + j 2 - h + fa 

a 4 = i 2 - i 3 + 2i 8 - 2i 9 + j 2 - j 3 + l 3 ~ h + h ~ h + k - fa + r 3 - r 4 + r 2 - r 3 

a 5 = r 1 -r 2 + i 1 -i 2 + k-l 7 + l 3 -k + 2i 7 - 2z 8 + r 4 - r 5 + ji - j 2 + fa - k 3 

a6 = h-k + i-h + h-h + Zie- ^7 + «"5 - r 6 + r 3 - r 4 + j - ji + k 2 -k 3 

a 7 = l 2 — l 3 + l 7 + h — i 2 + r 7 + r — r\ + j — ji + fa — fa + 2i 5 - 2i 6 

a 8 = I - h + / 4 - l 5 + i - it + r 6 - r 7 + r 2 - r 3 + 2z 4 - 2i 5 + fa - k 2 + ji - j 2 



and for the coherence of notations, (io,jo, ko, lo, ro) = (i,j,k,l,r) has to be under heard. Making use 
of the well-known formula 

8 

n\Si n 



dx r 



x 



x=0 



(A2) 



we end up with the form 



AiA 



1A 2 



7T 



-n\\n 2 \ 



i 



11 «e-l 

n e r^.i-v)! 

i,j,k,l,r=0 \ e=li e =0 



J E 



3 ie-l 



3 fee-l 

n e 



i 




7 r e -i 

n e 



i 



(r e _i-r e ) 



2«V 



2™ 



e=lr e =0 

l ,ni &b2,n 1^3,711 ^f)4,ni ^5 ,ri2 "&6,n2 ^bj ,n2^b%,ri2 



e=lfc e =0 

' x 2£ x r 

p 



(fee-l — fce)! 



(A3) 



*ll!i3^3^7^7! 

This shows clearly that a non vanishing purity should satisfy a set of constraint on different quantum 
numbers. These are 

61 — m = 2zn + i 3 - 14 - l 2 + I - r 2 + r + j 3 + k - fa - n\ = 

6 2 - m = 2«io - 2in + « 2 - i 3 + / 6 - r 7 + r 5 + j 3 + k 3 - m = 

63 - = 2« 9 - 2i w + i 3 - u ~ k + h + r 7 + r 4 - r 5 + j 2 - j 3 + fe 3 - m = 

6 4 - ni = i 2 - i 3 + 2z 8 - 2i 9 + j 2 — j 3 — Z 4 + Z 2 + k — fa — r 4 + r 2 — n x = 
b 5 -n 2 = l 5 -l 6 + i-i 1 + l 1 -l 2 + 2i 6 - 2i 7 + r 5 - r 6 + r 3 - r 4 + j - j x + fa - k 3 - n 2 = 
^6 — n 2 = l 2 — l 3 + l 7 + ii — i 2 + r7 + r — r\ + j — ji + fci — fa + 2i 5 - 2i 6 - n 2 = 
&7 - n 2 = n - r 2 + ii - i 2 + h ~ h + h ~ k + 2i 7 - 2i 8 + r 4 - r 5 + ji - j 2 + fa - k 3 - n 2 = 
bs — n 2 = l — h + k — l 5 + i — ii+r e — r 7 + r 2 — r 3 + 2i4 — 2i 5 + fa — k 2 + j\ - j 2 - n 2 = 0. 

(A4) 

We arrange the labels into two sets that we refer to them as the principals and secondary ones, 
respectively. The so-called secondary ones disappear upon summation of the 8 constraints and we get 



i + j + k + l + r = 2(ni + n 2 ) 



(A5) 



which is the constraint on the principal labels. The main result that emerges is that the purity is only 
depending on two parameters, such as 

P n , n Av,9)= . I "!,?' , a E C nin2 (i,j,fal,r)uVw k t l s r . (A6) 



sin 6» A y2cosh(2r;)+tan 2 |+cot 2 | 



2 i+i+fc+i+r=2(ni+ri2) 
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The most important future of our result is that the function C nin2 (i,j,k,l,r) can now be derived 
exactly for any n± and n 2 . This is 

/ 11 ie-l 1 \ / 3 je-l 1 \ / 3 fce-1 1 \ / 7 h-1 1 \ 

^e=li e =0 (*e-l - *e) ! / \e=lj e =0Ue-l ~ Jej! y y e =lfe e =0 ( fe e-l ~ fc e )! y \e=li e =0 ('e-1 ~ <e)" / 

(e=7r e _i ^ \ / 2~ i4 ( — I)* 2- * 8 (_]^ r_r i+ r 3 _r 5+r'6-r , 7 f-^\h—h+U—h+l&~h \ 

e 5 lr 5o(r e _i-r e )!j ( J ' (A?) 

Using the above constraints, we show that C niri2 can be reduced to the form 

4 ie-l\ / 3 je-l\ / 3 fee-l\ / 7 ie-l 1 \ / 7 



/ 4 i e -i \ / 3 J e -i \ / 3 fe e _i \ / 7 l e -! \ / 7 r e _i \ 

Cn in2 =nE nE nE nE n , n e t— rV 

\e=li e =0/ \e=lj e =0/ \e=lfe e =0/ \ e=ll e =0 \ l e-l ~ L e) ' / \e=lr e =0 l r e-l ~ r eJ ! / 

(A8) 



11 11 ^ 11 11 7T 

v e=li e =oy y e =lj e =0y ye=lfc e =0/ \e=ll e =0\ t e 

2~U f-j\h—ls+U—ls+l6—l7+r—ri+r3—r&+r6—r7+i2—ci—ca 



(i - h)\ (h - i 2 )l (i2 - kV- («3 - «4)!^W!c3!c4!c5!c6!c7!c 8 !c 9 !cio! 
where the involved parameters are fixed as 

ci = - [2ni - (i 2 - «3) - 0*2 - J3) - (A: - fei) - ('2 - h) - (h - k) - (r 2 - r 3 ) - (r 3 - r 4 ) 

(«3 - u) - (r-4 - r 5 ) - (j 2 - J3) - (k ~ h) ~ (k ~ k) ~r 7 - fc 3 ] 

C2 = - [ni - (i 2 - 13) - ( r 5 - r 6 ) - (r 6 - r 7 ) - / 6 - j 3 - A; 3 + (i 3 - u) - (I - h) 

(h - h) - (r - n) - (n - r 2 ) - j 3 - (A; - h)] 
c3 = , ni - (i 3 - n) - (r 4 - r 5 ) - (j 2 ~ J3) ~ (^4 ~ k) ~ (k ~ k) ~r 7 -k 3 ( ( 

C4 
C5 
C6 
C7 
C 8 
C9 
ClO 



2 

n\ - (i 2 - ^3) ~ (^5 ~ r 6 ) - (r 6 - r 7 ) - l e - j 3 - fc 3 

2 

rci ~ fe ~ u) ~ {I ~ h) ~ (h ~ h) ~ (r - n) - (n - r 2 ) - j 3 - (fc - fci) 

2 

«i - (i 2 - «3) - (J2 - J3) - (A; - fei) - (l 2 - h) ~ (h ~ k) - (r 2 - r 3 ) - (r 3 - r 4 ) 



n 2 


-(n 


- r 2 ) - (h 


- «2) - (^6 - "7] 
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I - (l 3 - i 4 ) 


- {r 4 


~r 5 ) 


-(Ji 
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~(k 2 
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*l) - («1 - fe) " 
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